Coupled dark energy: Towards a general description of the dynamics 
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In dark energy models of scalar-field coupled to a barotropic perfect fluid, the existence of 
cosmological scaling solutions restricts the Lagrangian of the field tp to p = Xg(Xe Xv ), where 
X = —g f *"d f i<pdvtp/2, A is a constant and g is an arbitrary function. We derive general evolution 
equations in an autonomous form for this Lagrangian and investigate the stability of fixed points for 
several different dark energy models-(i) ordinary (phantom) field, (ii) dilatonic ghost condensate, 
and (iii) (phantom) tachyon. We find the existence of scalar-field dominant fixed points (Q, v = 1) 
with an accelerated expansion in all models irrespective of the presence of the coupling Q between 
dark energy and dark matter. These fixed points are always classically stable for a phantom field, 
implying that the universe is eventually dominated by the energy density of a scalar field if phantom 
is responsible for dark energy. When the equation of state w v for the field <p is larger than —1, we 
find that scaling solutions are stable if the scalar-field dominant solution is unstable, and vice versa. 
Therefore in this case the final attractor is either a scaling solution with constant Q v satisfying 
< fitp < 1 or a scalar-field dominant solution with Q. v — 1. 

PACS numbers: 98.70.Vc 



I. INTRODUCTION 



Over the past few years, there have been enormous efforts in constructing models of dark energy-either motivated 
by particle physics or by phenomenological considerations (see Refs. [lj for review). The simplest explanation of dark 
energy is provided by cosmological constant, but the scenario is plagued by a severe fine tuning problem associated 
with its energy scale. This problem can be alleviated by considering a scalar field with a dynamically varying equation 
of state. In the recent years, a host of scalar field dark energy models have beenproposed, ranging from quintessence 
Q, k-essence 0, Born-Infeld scalars 0, phantoms [EEli ghost condensates 0,11] etc. 

In a viable dark energy scenario we require that the energy density of scalar field remains subdominant during 
radiation and matter dominant eras and that it becomes important only at late times to account for the current 
acceleration of universe. In this sense cosmological scaling solutions can be important building blocks in constructing 
the models of dark energy [j| 0, 0, 0] . The energy density of a scalar field ip decreases proportionally to that 
of a barotropic perfect fluid for scaling solutions. Steep exponential potentials give rise to scaling solutions for a 
minimally coupled scalar field in General Relativity allowing the dark energy density to mimic the background fluid 
during radiation or matter dominant era jl3l |. If the field potential V((p) becomes less steep at some moment of time, 
the universe exits from the scaling regime and enters the era of an accelerated expansion [l3l IhH ] . 

The quantity, A = —V V JV, which characterizes the slope of the potential, is constant for exponential potentials. In 
this case it is straightforward to investigate the stability of critical points in phase plane 9] . Even for general potentials 
a similar phase space analysis can be done by considering "instantaneous" critical points with a dynamically changing 
A 01 • Therefore we can understand the basic structure of the dynamics of dark energy by studying the fixed points 
corresponding to scaling solutions. Note that the potentials yielding scaling solutions are different depending upon 
the theories we adopt. For example we have V(ip) oc t^ -2 /^ 1-1 ) |12| for the universe characterised by a Friedmann 
equation: H 2 cx p n [the Randall- Sundrum (RS) braneworld and the RS Gauss-Bonnet braneworld correspond to n = 2 
and n = 2/3, respectively]. The scaling solution for tachyon corresponds to the inverse square potential: V(ip) oc ip~ 2 

US HUGH 

Typically dark energy models are based on scalar fields minimally coupled to gravity and do not implement the 
explicit coupling of the field to a background fluid. However there is no fundamental reason for this assumption 
in the absence of an underlying symmetry which would suppress the coupling. The possibility of a scalar field tp 
coupled to a matter and its cosmological consequences were originally pointed out in Refs. [18|. Amendola proposed a 
quintessence scenario coupled with dark matter [l9l | as an extension of nonminimal coupling theories |2(il ] (see Ref. pl| 
for an explicit coupling of a quintessence field to fermions or dark matter bosons). It is remarkable that the scaling 
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solutions in coupled quintessence models can lead to a late-time acceleration, while this is not possible in the absence 
of the coupling. Recently there have been attempts to study the dynamics of a phantom field coupled to dark matter 

mm 

In Refs. 0, [Q it was shown that the existence of scaling solutions for coupled dark energy restricts the form of 
the field Lagrangian to be p(X, <p) = X g(Xe Xv ), where X = —g^d^(pd u ip/2 and g is any function in terms of Xe Xv . 
This result is very general, since it was derived by starting from a general Lagrangian p(X, ip) which is an arbitrary 
function of X and ip. In fact this Lagrangian includes a wide variety of dark energy models such as quintessence, 
phantoms, dilatonic ghost condensates and Born-Infeld scalars. While the critical points corresponding to scaling 
solutions were derived in Ref . 0, ^3 1 this is not sufficient to understand the properties of all the fixed points in such 
models. In fact scaling solutions correspond to the fractional energy density of scalar fields satisfying < Q v < 1, 
but it is known that Q, v approaches 1 for several fixed points in the cases of an ordinary field ||| and a tachyon field 
[Tfil Hii Il7| when the coupling Q is absent between dark energy and dark matter. 

Our aim in this paper is to study the fixed points and their stabilities against perturbations for coupled dark energy 
models with the field Lagrangian: p(X,<p) — X g(Xe Xip ). This includes quintessence, dilatonic ghost condensate and 
tachyons, with potentials corresponding to scaling solutions. We shall also study the case of phantoms with a negative 
kinematic term in order to understand the difference from normal scalar fields. While phantoms are plagued by the 
problem of vacuum instability at the quantum level p^j , we would like to clarify the classical stability around critical 
points. We note that this quantum instability for phantoms is overcome in the dilatonic ghost condensate scenario 
provided that higher-order derivative terms stabilize the vacuum § . 

The rest of the paper is organised as follows. In Sec. II we briefly review the formalism of a general scalar 
field ip coupled to barotropic fluid and establish the autonomous form of evolution equations for the Lagrangian 
p(X,tp) — X g(Xe Xv ). In Sec. Ill, we apply our autonomous equations to the system with a standard (phantom) 
scalar field, obtain all the critical points and investigate their stabilities. Sec. IV and Sec. V are devoted to the 
detailed phase space analysis of dilatonic ghost condensate and (phantom) tachyon field, respectively. In Sec. VI we 
bring out some generic new features of coupled dark energy scenarios. 



II. SCALAR-FIELD MODEL 



Let us consider scalar-field models of dark energy with an energy density p and a pressure density p. The equation 
of state for dark energy is defined by w v = pj p. We shall study a general situation in which a field ip responsible for 
dark energy is coupled to a barotropic perfect fluid with an equation of state: w m = p m j ' p m . 

In the flat Friedmann-Robertson- Walker background with a scale factor a, the equations for p and p rn are 8] 

p + 3H(l + w v )p = -Qpmip, (1) 
p m + 3H(l + w m )p m = +Qp m (p , (2) 

where H = a/a is the Hubble rate with a dot being a derivative in terms of cosmic time t. The equation for the 
Hubble rate is 

H = -\[{l+w v )p+{l + w m )p m ], (3) 
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together with the constraint 



3H 2 = p + p m ■ (4) 



Here we used the unit 8irG = 1 (G is a gravitational constant). In Eqs. Q and © we introduced a coupling Q 
between dark energy and barotropic fluid by assuming the interaction given in Ref. [19(. In Refs. 22, 23j the authors 
adopted different forms of the coupling. In what follows we shall restrict our analysis to the case of positive constant 
Q > 0, but it is straightforward to extend the analysis to the case of negative Q. 

We define the fractional density of dark energy and barotropic fluid, Sl v = p/(3H 2 ) and f2 m = p m /(3H 2 ), with 
Cl v + Q, m = 1 by Eq. Scaling solutions are characterized by constant values of w v and fl v during the evolution. 
Then the existence of scaling solutions restricts the form of the scalar-field pressure density to be @, 

p = Xg(Xe^). (5) 



where X = —g^d^ipd u ip/2 and g is any function in terms of V = Xe Xv . Here A is defined by 

(6) 



w m - Q v {w m - w v ) 
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which is constant when scaling solutions exist. A is related with the slope of the scalar-field potential V(<j>). For 
example one has A oc ~V V /V for an ordinary field and A cx — V v /V 3 ^ 2 for a tachyon field 01 • Then the 
associated scalar-field potentials are given by V — V$e~ Xv for the ordinary field and an inverse power-law potential 
V = Vo<fi~ 2 for the tachyon. 

For the Lagrangian (JSJ the energy density for the field tp is 

P = 2*|| - P = X [g(Y) + 2Yg'(Y)\ , (7) 

where a prime denotes the derivative in terms of Y. Then Eq. Q can be rewritten as 

[g(Y) + Wg'(Y) + 2Y 2 g "(Y)\ ^ + ZH [g(Y) + Yg'(Y)] <p + XXY [3g'(Y) + 2Yg"(Y)] = -Qp m . (8) 
We introduce the following dimensionless quantities: 

y = /77„ • ( 9 ) 



V6H ' V3H 

From this definition, y is positive (we do not consider the case of negative H). 

Defining the number of e- folds as N = In a, we can cast the evolution equations in the following autonomous form: 



da; 
dN 



3 :(1 + w m ) + 1(1 w m )x 2 g(Y) - 3w m x 2 Yg'(Y) Y) " V)1 



2 y m ' 2 y m > * y ' m * y ' g(Y) + 5Yg'(Y) + 2Y 2 g"(Y) 
_ ^ \x 2 Y{Zg'(Y) + 2Yg"(Y)} + Q[l - x 2 {g(Y) + 2Yg'(Y)}} 

2[g(Y) + 5Yg>(Y)+2Y2g»(Y)} ' 1 } 

^ = -^Y^V + lv[ 1 + w ™ + - ^n)x 2 g(Y) - 2w m x 2 Yg'(Y)] , (11) 
= + Wm + (1 ~ Wm ^3( Y ) ~ 2w m x 2 Yg'(Y)] . (12) 



We also find 



Q v =x 2 \g(Y) + 2Yg'(Y)}, Wv = —— £^Q - - . (13) 
•p yv ;J ' 1/3 g(Y) + 2Yg'(Y) V ' 



It is also convenient to define the total effective equation of state: 



w cS ee ^ w ™ + (1 - ™ m )x 2 5 (F) - 2 V 2 y</(Y) . (14) 

P + Prn 



Combining Eq. l(T2l with Eq. ifHI) we obtain 

a 1 + 3ui e ff 

alP ~ 2 



(15) 



This means that the universe exhibits an accelerated expansion for w e Q < —1/3. 

It may be noted that the quantity Y is constant along a scaling solution, i.e., Y = Yq. However Y is not necessarily 
conserved for other fixed points for the system given by Eqs. I|10|) and (|llfl . Therefore one can not use the property 
Y = Yq in order to derive the fixed points except for scaling solutions. In subsequent sections we shall apply the 
evolution equations (|10fl and (|11|) to several different dark energy models- (i) ordinary (phantom) scalar field, (ii) 
dilatonic ghost condensate and (iii) (phantom) tachyon. 

By using Eqs. ©, (T3J|, and JT3J|, one can show that w e ff is written as 

w m X-Q 

WcS = ~qTY~ ' (16) 

whose form is independent of the function g(Y). For the pressureless fluid (w m = 0) we have w e e = —Q/(Q + A). 
Then we obtain w c g — > — 1 in the limit Q A. 

From Eq. l|13fl we find that the Q-dependent term in Eq. (|10fl drops out when Q, v approaches 1. Therefore, if the 
fixed point corresponding to Q v — 1 exists for Q = 0, the same fixed point should appear even in the presence of the 
coupling Q. This is a general feature of any scalar field system mentioned above and would necessarily manifest in 
all models we consider in the following sections. 
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III. ORDINARY (PHANTOM) SCALAR FIELD 

It is known that a canonical scalar field with an exponential potential 

p(X,<p) =eX-ce-^, (17) 

possesses scaling solutions. We note that e = +1 corresponds to a standard field and e = — 1 to aphantom. In fact 
this Lagrangian can be obtained by starting with a pressure density of the form: p = f(X) — V{(p) We obtain 

the Lagrangian (|17|l by choosing 

g(Y) = e-c/Y, (18) 

in Eq. 0. In what follows we shall study the case with A > without the loss of generality, since negative A 
corresponds to the change ip — ► — ip in Eq. I|17l) . 
For the choice itTHj) Eqs. lfTT)|) and (fTT|> reduce to 

^ = -3x + ^e\cy 2 + \x [(1 - w m )ex 2 + (1 + w m )(l cy 2 )} ^e(l - ex 2 - cy 2 ) , (19) 

^ - -^y\xy + \y [(1 - w m )ex 2 + (1 + w m )(l - cy 2 )] , (20) 

where Y is expressed through x and y, as F = x 2 jy 2 . We note that these coincide with those given in Ref. [j| for 
e = +1 and Q = 0. Eqs. JT3|) and (JT1} give 

2 2 

/"\ 2 2 / \ 2 / \2 / \ 

S2 V = ex +cy , w v = — — -, w e ff = w m + (1 - w m )ex - (1 + w m )cy . (21) 



A. Fixed points 

The fixed points can be obtained by setting dx/diV = and dy/dN = in Eqs. (|19fl and J3UJ. These are presented 
in Table I. 

• (i) Ordinary field (e = +1) 

The point (a) gives some fraction of the field energy density for Q ^ 0. However this does not provide an 
accelerated expansion, since the effective equation of state w c s is positive for < w m < 1. The points (bl) and 
(b2) correspond to kinetic dominant solutions with Q v = 1 and do not satisfy the condition w e f{ < —1/3. The 
point (c) is a scalar- field dominant solution (fl v = 1), which gives an acceleration of the universe for A 2 < 2. 
The point (d) corresponds to a cosmological scaling solution, which satisfies w v = w m for Q = 0. When 
the accelerated expansion occurs for Q > A(l + 3w m )/2. We note that the points (bl), (b2) and (c) exist 
irrespective of the presence of the coupling Q, since £l v = 1 in these cases. 

• (ii) Phantom field (e = —1) 

It is possible to have an accelerated expansion for the point (a) if the condition, Q 2 > (1 — w m )(l + 3w m )/2, is 
satisfied. However this case corresponds to an unphysical situation, i.e., £l v < for < w m < 1. The critical 
points (bl) and (b2) do not exist for the phantom field. Since w g = — 1 — A 2 /3 < —1 for the point (c), the 
universe accelerates independent of the values of A and Q. The point (d) leads to an accelerated expansion for 
Q > A(l + 3w m )/2, as is similar to the case of a normal field. 



B. Stability around fixed points 

We study the stability around the critical points given in Table I. Consider small perturbations u and v about the 
points (x c ,y c ), i.e., 

x = x c + u, y = y c + v- (22) 
Substituting into Eqs. (|10fl and (|ll|l . leads to the first-order differential equations for linear perturbations: 

£(:)-*(:)■ < 23 > 
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Name 


X 


y 


On? 1 




Weft 


(a) 







2Q 


1 


ir 1 2< ^ 


3e(l_ Wm ) 


3e{l — w m ) 


<"m T 3 e (l_,„ m ) 


(hi t 


1 

^r— 


n 
u 


1 

_L 


1 
r 


1 


(b2) 


•Si 





l 


1 


1 


(c) 


cX 




l 


1 ^ 3 


-1 + — 

1 ~ 3 


(d) 


•Se(i+w m ) 


r2Q(A+Q)+3e(l-< 1 ) 1 i/ 2 


Q(A+Q)+3e(l+™ m ) 


-Q(A+Q)+3€ju m (l+uj m ) 


™ m A-Q 


2(\+Q) 






Q(A+Q)+3e(l+w m ) 





TABLE I: The critical points for the ordinary (phantom) scalar field. The points (bl) and (b2) do not exist for the phantom 
field. 



where M. is a matrix that depends upon x c and y c . The elements of the matrix M. for the model (|17f) are 

9 3 

an = -3 + -ex 2 c {\ - w m ) + -(1 + w m )(l - cy 2 c ) + V6Qx c , (24) 
ai2 = \/6ecAy c - 3c(l + w m )x c y c + V6ecQy c , (25) 
a 21 = -^-\y c + 3ex c y c (l-w m ), (26) 

a-22 = -^Xx c -3c(l + w m )y^ + ^(l-w m )exl + ^(l + w m )(l-cy 2 ). (27) 

One can study the stability around the fixed points by considering the eigenvalues of the matrix M. . The eigenvalues 
are generally given by 



Mi, 2 



l r 



an + a 2 2 ± \J (an + a 22 ) 2 - 4(ouo 2 2 - ai2a 2 i] 



(28) 



We can evaluate /ii and /x 2 for each critical point: 
• Point (a): 

Mi = —77(1 - ™m) + -zj r 

2 e(l-w m ) 



1 



M2 



e(l - tu m ) 



Q(A + Q) + y(l-«&) 



Point (bl): 



V6. 



Mi = 3 —A, ^2 = 3(1 - w m ) + VoQ . 



• Point (b2): 



V6. 



Mi =3+— A, /i 2 = 3(1 - w m ) - V6Q . 



Point (c): 



Point (d): 



Mi = - 6 ) » M2 = eA(A + Q) - 3(1 + w m ) . 



Mi, 2 



3{A(1 - w m ) + 2Q} 
4(A + Q) 



1±W1 



8[3(1 + w m ) - eA(A + Q)] [3e(l - ta»,) + 2Q(A + Q)] 



3{A(1 - w m ) + 2QY 



One can easily verify that the above eigenvalues coincide with those given in Ref. for e = +1 and Q = 0. The 
stability around the fixed points can be generally classified in the following way: 
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• (i) Stable node: fi% < and /i2 < 0. 

• (ii) Unstable node: [i\ > and fi2 > 0. 

• (iii) Saddle point: fix < and \ii > (or /ii > and /i2 < 0). 

• (iv) Stable spiral: The determinant of the matrix Ai, i.e., T> = (an + CI22) 2 — 4(ana22 — 012021), is negative and 
the real parts of \i\ and 112 are negative. 

In what follows we shall analyze the stability around the fixed points for the ordinary field and the phantom. 

1. Ordinary field (e = +1) 

When the coupling Q is absent, the system with an ex pon ential potential was already investigated in Ref. 0. When 
Q =/= 0, the stability of fixed points was studied in Ref. [lJj when dust and radiation are present. We shall generally 
discuss the property of fixed points for a fluid with an equation of state: < w m < 1. 

• Point (a): 

In this case \X\ is negative if Q < y / 3/2(l — w rn ) and positive otherwise. Meanwhile /Z2 is positive for any value 
of Q and A (we recall that we are considering the case of positive Q and A). Therefore this point is a saddle for 
Q < ^/3/2(l — w m ) and an unstable node for Q > -\/3/2(l — w m ). We note that the condition < 1 gives 
Q < v/ (3/2)(l — w m ). Therefore the point (a) is a saddle for w m = under this condition. 

• Point (bl): 

While fj,2 is always positive, fii is negative if A > V6 and positive otherwise. Then the point (bl) is a saddle for 
A > \/6 and an unstable node for A < y/6. 

• Point (b2): 

Since [i\ is always positive and ^2 is negative for Q > (3/2) 1 / 2 (l — w m ) and positive otherwise, the point (c) is 
either saddle or an unstable node. 

• Point (c): 

The requirement of the existence of the point (c) gives A < \/6, which means that \i\ is always negative. The 
eigenvalue /12 is negative for A < (yf Q 2 + 12(1 + w m ) — Q)/2 and positive otherwise. Therefore the point (c) 
presents a stable node for A < {\/Q 2 + 12(1 + w m ) — Q)/2, whereas it is a saddle for {\/Q 2 + 12(1 + w m ) — 
Q)/2 < A < v/6. 

• Point (d): 

We first find that — 3{A(1 — w m ) + 2Q}/4(A + Q) < in the expression of \x\ and ^2- Secondly we obtain 
A(A + Q) > 3(1 + w m ) from the condition, fl v < 1. Then the point (d) corresponds to a stable node for 
3(1 + w m )/X — A < Q < and a stable spiral for Q > Q*, where Q* satisfies the following relation 

3 [A(l - w m ) + 2Q»] 2 = 8 [A(A + Q») - 3(1 + w m )\ [2Q»(A + Q„) + 3(1 - w 2 m )} . (29) 
For example we have Q* = 0.868 for A = 1.5 and w m = 0. 

The stability around the fixed points and the condition for an acceleration are summarized in Table II. The scaling 
solution (d) is always stable provided that fi v < 1, whereas the stability of the point (c) is dependent on the values 
of A and Q. It is important to note that the eigenvalue \i2 for the point (c) is positive when the condition for the 
existence of the point (d) is satisfied, i.e., A(A + Q) > 3(1 + w m ). Therefore the point (c) is unstable for the parameter 
range of Q and A in which the scaling solution (d) exists. 

In Fig. we plot the phase plane for A = 1.5, w m = and c — 1 with two different values of Q. In the phase space 
the allowed range corresponds to < x 2 + y 2 < 1. When A = 1.5 the point (c) is a saddle for Q > 0.5, and the point 
(d) is a stable node for 0.5 < Q < 0.868 and a stable spiral for Q > 0.868. The panel (A) in Fig. corresponds to 
the phase plane for Q = 0.6, in which case the point (d) is a stable node. We find that all trajectories approach the 
stable node (d), i.e., x c — 0.5832 and y c = 0.8825. In the panel (B) of Fig.^we plot the phase plane for Q = 4.0. It 
is clear that the critical point (d) [x c = 0.2227 and y c — 0.8825] is a stable spiral as estimated analytically. 
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Name 


Stability 


Acceleration 


Existence 


(a) 


Saddle point for Q < (3/2) 1/2 (l - w m ) 
Unstable node for Q > (3/2) 1/2 (l - w m ) 


No 


Q< (3/2) 1 / 2 (l-u- m ) 1 / 2 


(bl) 


Saddle point for A > \/6 
Unstable node for A < \/Q 


No 


All values 


(b2) 


Saddle point for Q > (3/2) 1/2 (l - w m ) 

Unstable node for (1 <T C3/9') 1 / 2 ('1 — in \ 


No 


All values 


(c) 


Saddle point for ([Q 2 + 12(1 + w m )] 1/2 - Q)/2 < A < ^6 
Stable node for A < ([Q 2 + 12(1 + w m )] 1/2 - Q)/2 


A < y/2 


A < \/6 


(d) 


Stable node for 3(1 + w m )/\ - A < Q < Q, 
Stable spiral for Q > Q t 


Q > A(l + 3w m )/2 


Q > 3(l+w m )/X- A 



TABLE II: The conditions for stability & acceleration & existence for an ordinary scalar field (e = +1). We consider the 
situation with positive values of Q and A. Here Q, is the solution of Eq. (1291 , 



(A) 




X 

FIG. 1: The phase plane for a standard scalar field corresponding to Q = 0.6 [panel (A)] and Q = 4.0 [panel (B)] for A = 1.5, 
w m ~ and c = 1. The late-time attractor corresponds to a stable node for Q = 0.6 and a stable spiral for Q = 4.0. The 
dotted curve is x 2 + y 2 = 1, which characterizes the border of the allowed region (fl v = 1). 

2. Phantom field (e = — 1) 

The fixed points (b) and (c) do not exist for the phantom field. 

• Point (a): 

In this case /n is always negative, whereas /j,2 can be either positive or negative depending the values of Q and 
A. Then this point is a saddle for Q(Q + A) < (3/2)(l - wf n ) and a stable node for Q(Q + A) > (3/2) (1 - w^). 
However, since VL V = — 2Q 2 /3(1 — w m ) < for < w m < 1, the fixed point (a) is not physically meaningful. 

• Point (c): 

Since both fix and [i-i are negative independent of the values of A and Q, the point (c) is a stable node. 

• Point (d): 

From the condition y 2 > 0, we require that 2Q(Q + A) > 3(1 — w^) for the existence of the critical point (d). 
Under this condition we find that [i\ < and fj,2 > 0. Therefore the point (d) corresponds to a saddle. 

The properties of the critical points are summarized in Table III. The scaling solution becomes always unstable for 
phantom. Therefore one can not construct a coupled dark energy scenario in which the present value of Q, v (~ 0.7) 
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Name 


Stability 


Acceleration 


Existence 


(a) 


Saddle point for Q(Q + A) < (3/2) (1 - u&) 
Stable node for Q(Q + A) > (3/2)(l - 


Q 2 > (1- w m )(l + 3w m )/2 


No if the condition < Sl v < 1 
is imposed 


(c) 


Stable node 


All values 


All values 


(d) 


Saddle 


Acceleration for Q > A(l + 3iu TO )/2 


Q(Q + \) > (3/2)(l-tujy 



TABLE III: The conditions for stability & acceleration & existence for a phantom scalar field (e = —1). We consider the 
situation with positive values of Q and A. 



1 . 5 




X 



FIG. 2: Phase plane for a phantom-type scalar field (e = —1) for Q = 3/4, A = 1, w m = and c = 1. The point A is a stable 
fixed point (c), whereas the point B is the fixed point (a) corresponding to a saddle. All trajectories approach the point A, 
which gives w e s = —4/3 and Q, v = 1. We also show the border of the allowed range (—1 < x 2 — y 2 < 0) as dotted curves. 

is a late-time attractor. This property is different from the case of an ordinary field in which scaling solutions can 
be stable fixed points. The only viable stable attractor for phantom is the fixed point (c), giving the dark energy 
dominated universe = 1) with an equation of state w v = —1 — A 2 /3 < — 1. 

In Figs. [21 and |3 we plot the phase plane for two different cases. Figure |21 corresponds to Q = 3/4, A = 1, w m = 
and c = 1. In this case the fixed point (a) is saddle, whereas the point (d) does not exist. It is clear from Fig. |3 
that the fixed point (c) is a global attractor. We also note that the allowed range in the phase plane corresponds to 
— l<x 2 — y 2 < 0, which comes from the condition < ft v < 1. The saddle point (a) exists outside of this region. 
Figure corresponds to Q — 3/2, A = 2, w m — and c = 1. In this case the fixed point (a) is a stable node, whereas 
there exists a saddle point (d). In Fig.|3|we find that the fixed points (a) and (d) are actually stable nodes, although 
the point (a) is not physically meaningful. Compared to Fig. the critical point (d) newly appears, but this is not 
a late-time attractor. It is worth mentioning that numerical results agree very well with our analytic estimation for 
the stability analysis about critical points. 

IV. DILATONIC GHOST CONDENSATE 

It was pointed out in Ref. |24| that at the quantum level a phantom field is plagued by a vacuum instability 
associated with the production of ghosts and photon pairs. We recall that the phantom field is characterised by a 
negative value of the quantity: px = dp/ dX. One can consider a scenario which avoids the problem of quantum 
instability by adding higher-order derivative terms that stabilize the vacuum so that px becomes positive 0- In 
the context of low-energy effective string theory one may consider the following Lagrangian that involves a dilatonic 
higher-order term [25| : 



p = -X + ce Xv X 2 , 



(30) 



9 




FIG. 3: Phase plane for a phantom-type scalar field (e = —1) corresponding to Q = 3/2, A = 2, w m = and c= 1. The point 
A is a stable fixed point (c), whereas the point B is the fixed point (a) corresponding to a stable node. The point C is the fixed 
point (d) corresponding to a saddle. In this case the trajectories approach either the point A or B, depending on the initial 
conditions of x and y. The border of the allowed range is the same as in Fig.|21 



where c and A are both positive. The application of this scenario to dark energy was done in Ref. |8(, but the stability 
of critical points was not studied. 

We obtain the pressure density (|3()[1 by choosing the function 



g(Y) = -l + cY. 



in Eq. ©. Then Eqs. (UJ and JTTJ yield 
dx 3 



dN 



dy_ 

dN 



x 1 1 + w m + (1 — w m )x (—1 + cY) — 2cw m x Y\ 



1 



1-QcY 
V6 



3(-l + 2cY)x + ^Xcx 2 Y + ^9. {l + x 2 (1 _ 3cF) } 



Axy + -y [l + w m + (1 - w m )x 2 (-l + cY) - 2cw m a; 2 Y] 



By Eqs. JTSJl and JTJJ we find 



= x 2 (— 1 + 3cY) , w v — 



-1 + cY 
-1 + 3cF 



, w cff = w m - (1 - w m )a; 2 + (1 - 3w m )cx 2 F . 



(31) 



(32) 
(33) 

(34) 



The stability of quantum fluctuations is ensured for px + 2Xpxx > and px > 0, which correspond to the 
condition cY > 1/2. In this case one has w v > —1 by Eq. I|34|l . which means that the presence of the term e Xv X 2 
leads to the stability of vacuum at the quantum level. 



A. Fixed points 



Setting dX/dN = and dY/dN = in Eqs. <|52" )l and <|53 | . we obtain 

2 _ 3cx 3 [A~ V6(Q + \)x] 2 = 3c(3w m - l)x 4 

V6x(\x - V&) + V6Q(l + x 2 ) ' V 3(l + w m )-3(l — w m )x 2 — v6Aa; 

Combining these relations, we get four critical points presented in Table IIVI The point (a) corresponds to x = 
and y = 0, in which case one has Y — > oo by Eq. I|35|l . Since we require the condition w m < — 1 in order to satisfy 
< f2,p < 1, this is not a physically meaningful solution. 
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Name 


X 


cY 








(a) 





oo 


3(^ m +l) 

3^ttj 1 


1/3 


-1 


(b) 


%/6A/ + (A) 
4 


1 , 

2 ' 16 


1 


-8+A^/_(A) 
8+3V/_(A) 


-8+A^/_(A) 
8+3A 3 /_(A) 


(c) 


\/6A/_(A) 
4 


1 | A a /+(A) 

2 i" 16 


1 


-8+A a / + (A) 
8+3A a / + (A) 


-8+A^/ + (A) 
8+3A^/ + (A) 


(d) 


y/B(l+w m ) 


3(l-ui^)-2Q(A+Q) 


3(l+iz> m )[l+w m -Q(A+Q)] 


3(l+™ m )iu m -Q(A+Q) 


w m \-Q 


2(A+Q) 


3(l-3ti> m )(l+w m ) 


(A+Q)*(l-3tB m ) 


3(l+» m )-3Q(A+(5) 


A+Q 



TABLE IV: The critical points for the ghost condensate model lIHOt . Here the functions /±(A) are defined by Eq. pijl . 



The points (b) and (c) correspond to the dark-energy dominated universe with fl v = 1. Therefore these points 
exist irrespective of the presence of the coupling Q. In Table ITvl the functions /±(A) are defined by 



/±(A) = 1± V /1 + 16/(3A 2 ). 



(36) 



Since one has w g 



■_p = (— l + cY)/(— l + 3cY) for the points (b) and (c), the condition for an accelerated expansion, 
w e fi < —1/3, gives cY < 2/3. We note that w e s < —1 for cY < 1/2 and w c s > —1 otherwise. The former case 
corresponds to the phantom equation of state, whereas the latter belongs to the case in which the system is stable 
at the quantum level. The parameter range of Y for the point (b) is 1/3 < cY < 1/2, which means that the field 
tp behaves as a phantom. In this case the universe exhibits an acceleration for any value of A and Q. The point (c) 
belongs to the parameter range given by 1/2 < cY < oo. The condition for an accelerated expansion corresponds to 
A 2 / + (A) < 8/3, which gives A < 0.817. In the limit A^Owe have cY — > 1/2, Q v — > 1 and w e n = w v — > — 1 for both 
points (b) and (c). The A = case is the original ghost condensate scenario proposed in Ref. 0, i.e., p = —X + X 2 , 
in which case one has an equation of state of cosmological constant (w v = — 1). 

The point (d) corresponds to a scaling solution. Since the effective equation of state is given by u> c ff = (w m \ — 
Q)/(A + Q), the universe accelerates for Q > A(l + 3w m )/2. The requirement of the conditions < ft v < 1 and 
Y = x 2 /y 2 > places constraints on the values of Q and A. For the non-relativistic dark matter (w m — 0) with 
positive coupling (Q > 0), we have the following constraint 



\/9Q 2 + 12 



< A< - - Q. 



(37) 



This implies that Q needs to be smaller than 1 for positive A. 



B. Stability around fixed points 

The elements of the matrix M. for perturbations u and v are given by 

-3 + V6x c Q + 18cY c + 6%/6c(A - Q)x c Y c 



an 



an 



a-22 



[I + w m - 3(1 - w m )x c + 5c(l - 3w m )x c Y c \ + 



12cK 



x c (-l + 6cY c y 



3(-l + 2cY c )x c 



3V6 



\cx C Y C 



\/6Q 



{1- 



1 - 6cY c 
l(-l + 3cY c )} 



oia = -3c(l - 3w m )x 2 c Y c 3 / 2 - 



V6c 



1 - 6cY c 



r c 3/2 {2V6 + 3x c (A-Q)} 



12cr r 



y c (-l + 6cY C/ 



3(-l + 2cY c )x c + ^Xcx 2 c Y c + ^ {1 - x\{-\ + 3cY c )} 



>/6, 



3y c [-(1 - w m )x c + 2c(l - 3w m )x c Y c ] Y^ Vc 

3 

- [l + w m - (1 - w m )x 2 c - c(l - 3w m )x 2 Y c 



V6 



(38) 



(39) 

(40) 
(41) 



Since the expression of the matrix elements is rather complicated, the eigenvalues of A4 are not simply written unlike 
the case of Sec. III. However we can numerically evaluate /ii and fi2 and investigate the stability of fixed points. 
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Name 


Stability 


Acceleration 


Existence 


(a) 


Unstable 


All values 


w m < —1 


(b) 


Stable spiral or stable node 


All values 


All values 


(c) 


Stable node or saddle 


A < 0.817 


All values 


(d) 


Stable node or stable spiral or saddle 


Q > \(l + 3w m )/2 


[(9Q 2 + 12) 1/2 - 5Q]/2 < A < 1/Q - Q for w m = 



TABLE V: The conditions for stability & acceleration & existence for the ghost condensate model 1301 . 



• Point (a): 

In this case the component an diverges, which means that this point is unstable in addition to the fact fi v does 
not belong to the range < Q v < 1 for plausible values of w m . 

• Point (b): 

We numerically evaluate the eigenvalues /xi and /i2 for w m = 0, c = 1 and find that the point (b) is either a 
stable spiral or a stable node. When Q < 10 the determinant of the matrix M. is negative with negative real 
parts of /ii and /j.2, which means that the point (b) is a stable spiral. In the case of Q > 10, the determinant 
is positive with negative /xi and /i2 if A is smaller than a value A*(<3), thereby corresponding to a stable node. 
Here A»(Q) depends on the value Q. When A > \*(Q) and Q > 10, the point (b) is a stable spiral. 

• Points (c) & (d): 

It would be convenient to discuss the stability of these critical points together as there exists an interesting 
relation between them. We shall consider the case of w m = and c = 1 . 

For the point (c) we have numerically found that \ii is always negative irrespective of the values of Q and A. 
Our analysis shows that there exists a critical value A*(Q) such that fix is negative for A < A*(Q) and becomes 
positive for A > A*(Q). The critical value of A can be computed by demanding (ana22 — 012021) = 0, which 
leads to 

1 



A*(Q) = - [V9Q2 + 12 - 5Qj . (42) 

We conclude that in the region specified by < A < A* (Q) , the critical point (c) is a stable node which becomes 
a saddle as we move out of this region [A > A*(Q)]. 

For the point (d) the second eigenvalue fj,2 is negative or Re (/i 2 ) < for all values of Q and A. Meanwhile the 
first eigenvalue exhibits an interesting behavior. We recall that the allowed domain for the existence of the point 
(d) lies outside the region of stability for the point (c), see Eq. QH7JI. If we extend it to the region A < A*(Q), we 
find that point (d) is a saddle in this domain (/ii > and (12 < 0). The critical value of A, at which fi\ for the 
point (d) vanishes, exactly coincides with A*(Q) given by Eq. I|42l) . As we move out of the domain of stability 
for the point (c), the critical point (d) becomes a stable node as \i\ < in this case. We numerically find there 
exists a second critical value A*»(Q) (> A*(Q)) at which the determinant T> of the system vanishes such that 
/ii < for A*(Q) < A < A**(Q) and Re(/Lti) < for A > A**(Q). To have an idea of orders of magnitudes, let 
us quote some numerical values of A*(Q) and A**(<2), for instance A**(Q) = 2.06, 1.15; A*(Q) = 1.26, 0.64 for 
Q = 0.2,0.5 respectively. The stability of (c) & (d) can be briefly summarised as follows. The point (c) is a 
stable node whereas the point (d) is a saddle for < A < A*(Q). The point (c) becomes a saddle for A > A»(Q). 
In the region characterised by A»(Q) < A < A**(Q), the critical point (d) is a stable node but a stable spiral for 
A > A*»(<3). 

We summarize the property of fixed points in Table IVI Although the point (b) is always stable at the classical level, 
this corresponds to a phantom equation of state (w v < — 1). Therefore this is plagued by the instability of vacuum at 
the quantum level, whereas the point (c) is free from such a quantum instability. The scaling solution (d) also gives 
rise to an equation of state w v > — 1 as can be checked by the expression of w v in Table IIVI When the point (c) is 
stable we find that the point (d) is unstable, and vice versa. Therefore the final viable attractor is described by a 
dark energy dominant universe with Q v = 1 [case (c)] or by a scaling solution with < f2 v — const < 1 [case (d)]. 

In Fig. 0] we plot the phase plane for Q — 0.5 and A = 0.4. In this case the point (c) is a saddle, whereas the point 
(d) is a stable node. By using the condition < ft v < 1 in Eq. (|34|l . we find that x and y are constrained to be in the 
range: 3cx 4 /(l + x 2 ) < y 2 < 3cx 2 . We also obtain the condition, y 2 < 2cx 2 , if the stability of quantum fluctuations 
is taken into account 0. When x and y are initially smaller than of order unity with positive x, the trajectories tend 
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FIG. 4: Phase plane analysis in the dilatonic ghost condensate case for Q = 0.5 and A = 0.8 with w m = and c = 1. In this 
case there exist two stable points (b) (x,y) = ( — 1.98, 3.07) and (d) (x,y) = (0.94, 1.25). The curves (i), (ii) and (iii) correspond 
to y 2 = 3cx 4 /(l + x 2 ), y 2 — 3cx 2 and y 2 — 6cx 2 , respectively. 



to approach the line y = \/6cx on which the speed of sound, c s = y/px/px, diverges js|. Therefore these cases are 
physically unrealistic. Meanwhile when initial conditions of x and y are not much smaller than of order unity, the 
solutions approach the stable point (d) provided that x is positive (see Fig.0J. When x and y are much smaller than 
1 during matter dominant era, it is difficult to reach the critical point (d) for constant Q. If the coupling Q rapidly 
grows during the transition to scalar field dominated era, it is possible to approach the scaling solution (d) [8| . 

If the initial value of a; is negative, the trajectories approach the stable point (b). In this case we numerically found 
that the solutions do not cross the lines y = y/&c\x\ even for the initial values of |a;| much smaller than unity. The 
final attractor point (b) corresponds to the phantom-dominant universe with Q v = 1. 



V. TACHYON AND PHANTOM TACHYON 

The Lagrangian of a tachyon field is given, in general, by [22| 

p=-V(<t>)y/l-e^, (43) 

where e = +1. Here we allowed the possibility of a phantom tachyon (e = —1) (2a. We note that there are many 
works in which the dynamics of tachyon is investigated in a cosmological context |29j . The potential corresponding 
to scaling solutions is the inverse power-law type, i.e., V(<fi) oc <p~ 2 . If we choose the function g(Y) as 

g(Y) = -cVl - 2eY/Y , (44) 

the Lagrangian © yields p = —ce~ Xip y/l — 2eY. With a suitable field redefinition given by <f) = (2/A)e A<p / 2 , we obtain 
the tachyon Lagrangian l|43|l with an inverse square potential: V{4>) = 4c/(A 2 (/> 2 ). Note that we are considering 
positive c. 

We shall introduce a new variable x which is defined by i 2 = cj) 2 /2 = Y. Then x and i are now related as x — xy. 
Noting that 

C6 

g(Y) + 5Yg'(Y) + 2Y 2 g"(Y) = (45) 
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Name 


X 


2/ 




w v 


Wefi 


(a) 











-1 






i 


u 


1 
1 


u 


u 


(b2) 


1 

s/2i 





1 








(c) 


^Vc 
\/6e 


Vc 


1 


3e X 


3e ± 


(dl) 


(7/2,) 1/2 


A 


3ce7 
A^Vl-7 


Wm 


Wm 


(d2) 


-(7/2e) 1/2 


A 


3ce7 




W m 



TABLE VI: The critical points for the tachyon model (1431 with Q = 0. Here j/ c is defined by Eq. 150H . The points (b) and (d) 
do not exist for the phantom tachyon (e < 0). 



we obtain the following equations by using Eqs. (|10fl and i|ll|) : 



where 

Eqs. (O and Jll|) give 

bit in 



;, ( _^ + V^ (A/r -^_ cy2) 



— - ~ — \xy +- 



cy 



VI - 2ei 2 



2e 2cey 

,2 



cy 



VI - 2ei 2 



= (7- 2ei^) 



7 = 1 + io„ 



w„ = 2ex 2 - 1 , w cf f = w, n 



cy 2 {l + w m - 2ex 2 ) 
VI - 2ei 2 



(46) 
(47) 

(48) 

(49) 



A. Fixed points 

When we discuss fixed points in the tachyon model, it may be convenient to distinguish between two cases: 
and Q ^ 0. We note that the fixed points were derived in Ref. 

HG3 for Q = and e = +1. 



1. Q = 

We summarize the fixed points in Table IVll for Q — 0. The point (a) is a fluid-dominant solution (il m — > 1) with an 
effective equation of state, w e s = w m . Then the accelerated expansion does not occur unless w m is less than —1/3. 
The points (bl) and (b2) are the kinematic dominant solution whose effective equation of state corresponds to a dust. 
This does not exist in the case of phantom tachyon. The point (c) is a scalar-field dominant solution that gives an 
accelerated expansion at late times for \ 2 y 2 /(3e) < 2/3, where y c is defined by 



-A 2 ± VA 4 + 36c 2 



6ec 2 



"I 1/2 



(50) 



This condition translates into A 2 < 2V3c for e = +1. In Eq. (|50() the plus sign corresponds to e > 0, whereas the 
minus sign to e < 0. Note that phantom tachyon gives an effective equation of state w e e that is smaller than — 1. 
The points (dl) and (d2) correspond to scaling solutions in which the energy density of the scalar field decreases 



proportionally to that of the perfect fluid (u> ¥ 



j). The existence of this solution requires the condition w m < 



as can be seen in the expression of We note that the fixed points (dl) and (d2) do not exist for phantom tachyon 
unless the background fluid behaves as phantom (w m < — 1). 
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2. Q^O 



Let us next discuss the case with Q ^ 0. The critical points are summarized in Table rVTll We first note that the 
fixed point (x, y) = (0, 0) disappears in the presence of the coupling Q. The points (b) and (c) appear as is similar to 
the case Q — 0. While (bl) and (b2) do not exist for phantom tachyon, the point (c) exists both for e > and e < 0. 

The point (d) [(x,, j/j)] corresponds to scaling solutions, whose numbers of solutions depend on the value of w m . yi 
are related with Xi through the relation y^ — 3 , f/[^/6(X + Q)xj\. Here Xi satisfy the following equation 



7 - 2ct 2 2Q(A + Q) 
xfy/1 - 2exi ~ 3(l + w m ) C - 

We note that only positive Xi are allowed for A > 0, since j/j is positive definite. 
= yj\ — 2exf, we find 



(51) 



If we introduce the quantity 



f(0 = 



e+W m = Q{Q + A) 

£(l-£ 2 ) 3ce(l + Wm ) 



(52) 



The behavior of the function /(£) is different depending on the value of w m . We note that A is related with the slope 
of the potential as A = —V^jV I 2 [l^, which means that A > for the inverse power-law potential: V{4>) = M 2 <f>~ 2 . 
Then the r.h.s. of Eq. 1)52(1 is positive for e > and negative otherwise. The allowed range of £ is < £ < 1 for e > 
and £ > 1 for e < 0. We can classify the situation as follows. 



• w m > 0: 

The function /(£) goes to infinity for £ 



-> + and g 

Here £m is defined by = [-(1 + 3w m ) + y/(l + w„ 
for Eq. (52J provided that Q(Q + A)/(3ce(l 
for £ ^ 1 + and /(£) -> for £ - 



1 — 0. It has a minimum at £ = £a/ with < £zvf < 1. 



)(1 + 9w m )]/2. Then for e > 0, there exist two solutions 
w m )) > /(^m)' The function /(£) has a dependence /(£) — > — oo 
oo. Then for e < 0, there exists one scaling solution for Eq. (52} < 



0: 



In this case we can analytically derive the solution for Eq. 152f) . The function /(£) is zero at £ = and 
monotonically increases toward +oo as £ — * 1 — 0. Then we have one solution for Eq. (|52|1 if e > 0. The function 
/(£) has a dependence /(£) — > — oo for £ — > 1 + and /(£) — > for £ — > +oo. This again shows the existence of 
one solution for Eq. I|52l) if e < 0. The solutions for Eq. (|52ll are given by £ = [— 1 + vt 
[1 



and £ = [1 + Vl + 4 J 4 2 ]/2A for e = -1, where A 
points together with the equation of state w^: 



Q(Q + X)/(3c(l + w, 



4A 2 ]/2A for e = +1 
t )). Then we obtain the following fixed 



V^/c(-3c + ^9c 2 + 4Q 2 (A + Q) 2 ) 

2Q(A + Q) 
V6Q 

V3y/c(-3c + y/9c 2 + 4g 2 (A + Q) 2 ) 
9c 2 



4Q 2 (Q + A) 2 



1 , 4Q 2 (Q + A) 2 
9c 2 



(53) 
(54) 

(55) 



for e = +1 and 



V3^/c(3c+ ^9c 2 + 4Q 2 (A + g) 2 ) 
2Q(A + Q) 
V6Q 



V3Vc(3 c + ^9c 2 + 4Q 2 (A + Q) 5 



9r- 



4Q 2 (Q + A) 2 



i | 4Q 2 (Q + A) 2 | i 



9c 2 



(56) 
(57) 

(58) 



for e = — 1. 
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Name 


X 


y 




w v 


Weff 


(bl) 


1 

\/2e 





1 








(b2) 


1 

•J2i 





1 








(c) 


\Vc 
V6e 


Vc 


1 


^'2 '2 

x vi i 
3e L 


V* - 1 


(d) 


X% 




3c 7 2 


2tx\ - 1 


ey jf (l-ftu m — 2e«jf ) 


V6{\ + Q)ii 


2(A+Q) 2 i 2 A /l-2ei 2 


l^m 1 tt 



TABLE VII: The critical points for the tachyon model with Q / 0. Here y c is defined by Eq. a;, are the solutions of 

Eq. 151H . The critical points (bl) and (b2) do not exist for the phantom tachyon (e < 0). The numbers of the point (d) depend 
on the value of w m . 



For ordinary tachyon one has x c — ► l/\/2, y c — ► -\/3/(A + Q), w v — > as Q — > and ai c — * 0, y c — * l/v 7 ^, 
u> v — + — 1 as Q — > +oo. There is a critical value Q*{\) which gives the border of acceleration and deceleration, 
i.e., the accelerated expansion occurs for Q > Q*(\). For phantom tachyon we obtain x c — ► oo, y c — ► 0, 
Wtp — > — oo as Q — > and x c — ► 0, y c — * l/y/c, w v — ► — 1 as Q — * +oo. Thus the presence of the coupling Q can 
lead to an accelerated expansion. 

w m < 0: 



The function /(£) is zero at £ = \f—w m . It monotonically increases toward +oo as £ — » 1 — in the region 
v 7 — w TO < £ < 1. This means that we have one solution for Eq. fB^f if e > 0. The function /(£) has a dependence 
/(£) — > — oo for £ — * 1 + and /(£) — * for £ — ► +oo, which shows the existence of one solution for Eq. (15211 if 
e < 0. Note, however, that the w m < case is not realistic. 



B. Stability 



The components of the matrix M. are 



an 



-3 + 18 ex 2 c - 2V6(A + Q)x c y c + 3 ^ Qxc ^1 - 2ei* , 



oia = (l-2e^) 



\/6(A + Q) V6Q 
2e 2ce2/2 



«21 = 



022 



"~2~ Vc 



(l- 2e^)3/2 yr - 2^f' 



/- 3 9cy c 2 ( 7 -26xg) 

-vdA.t c w c + —7 — 

" 2 2^1 - 2eJ2 



(59) 
(60) 
(61) 
(62) 



Hereafter we shall discuss the stability of fixed points for an ordinary tachyon (e = +1) and for a phantom tachyon 
(e = —1) by evaluating eigenvalues of the matrix M.. 



1. Ordinary tachyon (e = +1) 

The stability of fixed points is summarized in Table IVlTTl 
• Point (a): 

This point exists only for Q = 0. The eigenvalues are 

Mi = 3 7 /2 , [i 2 = -3 . (63) 

Therefore the point (a) is a saddle for 7 > and a stable node for 7 < 0. Therefore this point is not stable for 
an ordinary fluid satisfying 7 > 1. 
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• Points (bl) and (b2): 
Since the eigenvalues are 

Mi = 6 , /i 2 = 9/2 - 3 7 , (64) 
the points (bl) and (b2) are unstable nodes for 7 < 3/2 and saddle points for 7 > 3/2. 

• Point (c): 

The eigenvalues are 

Mi = -3 7 + A(A 6 ^ 2 Q) (Va 4 + 36c 2 - A 2 ) , ^ = -3+^(^ + 36^ -A 2 ) . (65) 
The range of M2 is —3 < M2 < —3/2. We also find that Mi < if 

7>7^^Q)(v / A T T^-A 2 ) , (66) 

and mi > if 7 < 7s . Therefore the point (c) is a stable node for 7 > 7s and a saddle point for 7 < 7s . 
Point (d): 

When Q = we obtain 7 < 7^ = A 2 [VA 4 + 36c 2 — A 2 ]/(18c 2 ) from the condition £l v < 1. The eigenvalues for 
Q = are 



Mi,2 = i 



48c 7 2 



7 - 2 ± ^/ 17 7 2 - 20 7 + 4 + — !- y/1 -7 



(67) 



which arc both negative for 7 < 7s . Therefore the point (d) is a stable node for Q = 0. 

This situation changes if we account for the coupling Q. In what follows we shall consider the case of a non 
relativistic dark matter (w m — 0). The analytic expressions for the eigenvalues are rather cumbersome for 
general Q, but they take simple forms in the large coupling limit (Q — > +00): 



Mi, 2 ^-3±y-^. (68) 

This demonstrates that the critical point (d) is a stable spiral for large Q. In fact we numerically confirmed 
that the determinant of the matrix M. changes from positive to negative when Q becomes larger than a critical 
value <3i(A). This critical value depends on A, e.g., Qi(X = 0.1) = 37.5 and <5i(A = 1.0) = 3.28 for c = 1. When 
Q > Qi(X) numerical calculations show that the point (d) is a stable spiral for any A. 

When Q < Qi(X) we find that both mi and M2 are negative when Q is larger than a critical value (52(A). 
Meanwhile Mi > and M2 < for Q < (32(A). Here Q2(A) can be analytically derived as 

A VA 4 + 36c 2 

Q 2 (A) = -- + ^-^ , (69) 

which corresponds to the eigenvalue: mi — 0- For example we have Q2(A = 0.1) = 29.95 and Q2(X = 1.0) = 2.54 
for c = 1. From the above argument the fixed point (d) is a saddle for Q < (52(A), a stable node for (52(A) < 
<5 < <5i(A) and a stable spiral for Q > Qi(A). 

In the case of w m — 0, the stability condition (|66|) for the point (c) corresponds to 



The r.h.s. completely coincides with Q 2 (A). This means that the critical point (c) presents a stable node in the region 
where (d) is a saddle. When Q > (52(A) the point (d) is stable, whereas (c) is a saddle. Therefore one can not realize 
the situation in which both (c) and (d) are stable. 

In Fig. [S] we plot the evolution of VL V and fi m for Q = 3.0 and A = 2.18. Since <5i(A) = 1.02 and (52(A) = 0.67 
in this case, the fixed point (d) is a stable spiral whereas the point (c) is a saddle. In fact the solutions approach 
the point (d) with oscillations as is clearly seen in Fig. |SJ The attractor corresponds to a scaling solution that gives 
tt v = 0.7 and Q m = 0.3. 
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i. 1 QjlllC 


Stahilitv 


A ppplpratinn 


1 J A 1 ii LdlVjC 


(a) 


Saddle point for 7 > 
Stable node for 7 < 


7 <2/3 




(bl), (b2) 


Unstable node for 7 < 3/2 
Saddle point for 7 > 3/2 


JNo 


All values 


(c) 


Stable node for 7 > 7s 
Saddle point for 7 < 7s 


X 2 < 2\/3c 


All values 


(d) [Q = 0] 


Stable node or stable spiral 


7<2/3 


< 7 < 1 


(d) [Q + 0] 


Stable node or stable spiral or saddle 


Q > Q,(A) 


x 2 (l n x 2 ) 1/2 ^ ^ 



TABLE VIII: The conditions for stability & acceleration & existence for an ordinary tachyon field (e = +1). 7s is defined by 
Eq. 16611 . Q*(A) depends on A. 

1 .0 



0.80 



0.60 



0.40 



0.20 



0.0 

II 2 4 6 8 1 

N 

FIG. 5: The evolution of Sl v and il m for an ordinary tachyon (e = +1) with Q = 3.0 and A = 2.18. The stable attractor in 
this case is the fixed point (d), giving Q.^, — 0.7 and Q m = 0.3 asymptotically. 

2. Phantom tachyon (e = — 1) 

For a phantom tachyon the stability of fixed points exhibits a number of differences compared to an ordinary 
tachyon. 

• Point (a): 

The property of this fixed point is completely the same as in the case of an ordinary tachyon. 

• Point (c): 

In this case the calculation of the eigenvalues of the matrix M is more involved compared to the point (c) for 
e = +1, but we numerically find that fi± and are both negative for any value of Q and A. Then the point (c) 
is a stable node. 

• Point (d): 

In the case of w m = 0, the critical points are given by Eqs. I|56l) and (|57|l . In the large coupling limit Q — > 00, 
the eigenvalues are approximately given by 

m,2 - -3± W— , (71) 

which means that fix is positive and ^2 is negative. Numerically we find that jix > and /12 < for any value 
of Q and A. Then the fixed point (d) is always saddle. 
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Name 


Stability 


Acceleration 


Existence 


(a) 


Saddle point for 7 > 
Stable node for 7 < 


7<2/3 


Q = 


(c) 


Stable node 


All values 


All values 


(d) 


Saddle point 


All values 


Q ± and x\{\ + > 



TABLE IX: The conditions for stability & acceleration & existence for a phantom tachyon field (e = — 1). x% are the solutions 
of Eq. 

From the above argument scaling solutions do not give a viable late-time attractor for the phantom. This property 
is similar to the ordinary phantom field discussed in Sec. III. Thus the solutions approach the phantom dominant 
universe (f^ = 1) even in the presence of the coupling Q. 

VI. SUMMARY 

In this paper we studied coupled dark energy scenarios in the presence of a scalar field tp coupled to a barotropic 
perfect fluid. The condition for the existence of scaling solutions restricts the form of the field Lagrangian to be 
p = Xg(Xe Xip ), where X = —g tlv d tl ipd u tp/2, A is a constant and g is an arbitrary function 0. This Lagrangian 
includes a wide variety of scalar-field dark energy models such as quintessence, diatonic ghost condensate, tachyon 
and k-essence |12| . Our main aim was to investigate, in detail, the properties of critical points which play crucially 
important roles when we construct dark energy models interacting with dark matter. 

We first derived differential equations (|10|) and (jllfl in an autonomous form for the general Lagrangian p = 
Xg{Xe Xv ) by introducing dimensionless quantities x and y defined in Eq. ©. These equations can be used for 
any type of scalar-field dark energy models which possess cosmological scaling solutions. We note that the quantity 
A is typically related with the slope of the scalar-field potential V(tp), e.g., A cx —V v /V for an ordinary field [llj and 
A cx -V v /V 3 / 2 for a tachyon field [13 ■ Scaling solutions are characterized by constant A, thus corresponding to an 
exponential potential V = Voe~ Xv for an ordinary field and an inverse power-law potential V = Voip~ 2 for tachyon. 
Even for general potentials one can perform a phase-space analysis by considering "instantaneous" critical points with 
a dynamically changing A [ill Il7| . Thus the investigation based upon constant A contains a fundamental structure 
of critical points in scalar-field dark energy models. 

We applied our autonomous equations to several different dark energy models-(i) ordinary scalar field (including 
phantom), (ii) dilatonic ghost condensate, and (iii) tachyon (including phantom). In all cases we found critical points 
corresponding to a scalar-field dominant solution (Q v — 1) with an equation of state: w v — > —1 as A — > 0. These 
points exist irrespective of the presence of the coupling Q. This can be understood by the fact that the Q-dependent 
term in Eq. (|10(l vanishes for tt v — > 1. In the case where w v > — 1, these solutions are either stable nodes or saddle 
points depending on the values of A and Q, see (d) in Table II, (c) in Table V, and (c) in Table VIII. We note that 
the condition for an accelerated expansion requires that A is smaller than a critical value A, i.e., A < v2 for the 
ordinary field, A < 0.817 for the dilatonic ghost condensate, and A 2 < 2-\/3c for the tachyon. The current universe 
can approach this scalar-field dominated fixed point with an accelerated expansion provided that this point is a stable 
node and A is smaller than A. In the case of a phantom field (w v < — 1), the Q-independent solutions explained above 
are found to be always stable at the classical level, see (d) in Table III, (b) in Table V, and (c) in Table IX. Thus the 
solutions tend to approach these fixed points irrespective of the values of A and Q. Nevertheless we need to keep in 
mind that this classical stability may not be ensured at the quantum level because of the vacuum instability under 
the production of ghosts and photon pairs H, . 

In the presence of the coupling Q there exist viable scaling solutions that provide an accelerated expansion, while 
it is not possible for Q = unless w m is less than —1/3. If w v > — 1, the scaling solution is a stable node or a stable 
spiral for an ordinary scalar field under the condition Q, v < 1, see (e) in Table II. In the cases of dilatonic ghost 
condensate and tachyon, the scaling solution can be a stable node or a stable spiral or a saddle depending on the 
values of A and Q, see (d) in Table V and (d) in Table VIII. The accelerated expansion occurs when the coupling Q 
is larger than a value Q*(A), e.g., Q > A(l + 3w m )/2 for the ordinary field and the dilatonic ghost condensate. When 
w v > — 1 we find that the scaling solution is stable if the scalar-field dominated fixed point — 1) is unstable, 
and vice versa. This property holds in all models considered in this paper. Therefore the final attractor is either the 
scaling solution with constant VL V satisfying < Cl v < 1 or the scalar-field dominant solution with Q v = 1. 

For the ordinary phantom and the phantom tachyon we found that scaling solutions always correspond to saddle 
points, see (e) in Table III and (d) in Table IX. Therefore they can not be late-time attractors unlike the case of 
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w v > — 1. The situation is similar in the dilatonic ghost condensate model as well, since the solutions do not reach a 
scaling solution for initial values of \x\ and y much smaller than 1. The only viable stable attractor is the scalar-field 
dominant fixed point corresponding to f2„ = 1 and < —1. Therefore the universe finally approaches a state 
dominated by the phantom field even in the presence of the coupling Q. This tells us how phantom is strong to over 
dominate the universe! 

Since our paper provides a general formalism applicable to a variety of scalar fields coupled to dark matter, we hope 
that it would be useful for the concrete model building of dark energy. Recently Amendola et al. [3(J carried out a 
likelihood analysis using the dataset of Type la supernovae for the model i|17fl arid placed constraints on the values 
of w e g and They found that the coupled dark energy scenario is compatible with the observational data only if 
the equation of state satisfies w c g > — 1. It would be quite interesting to extend the analysis to general coupled dark 
energy models presented in this paper. 
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